Unusual scaling for pulsed laser deposition 
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We demonstrate that a simple model for pulsed laser deposition exhibits an unusual type of scaling 
behavior for the island density in the submonolayer regime. This quantity is studied as function of 
pulse intensity and deposition time. We find a data collapse for the ratios of the logarithms of these 
quantities, whereas conventional scaling as observed in molecular beam epitaxy involves ratios of 
powers. 
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Many systems in equilibrium and nonequilibrium sta- 
tistical physics exhibit power-law scaling. This means 
that a system with an observable, M, depending for ex- 
ample on two parameters, z\ and z 2 , looks the same, if 
the units of M , z\ and z 2 are rescaled by certain factors 
which are related to each other by power-laws. Such a 
scaling transformation can be written as 

zi -> Azi, z 2 -> A z 2 , M -> A a M , (1) 

where A is a scaling parameter and a, (3 are certain ex- 
ponents. Eq. (|l|) implies the scaling form 

M(z 1 ,z 2 )=z?f(z 2 /zP), (2) 

where / is a scaling function depending on a scale- 
invariant argument. This type of scaling can be observed 
in a vast variety of applications, including equilibrium 
critical phenomena |l|,0), growth processes &-@], driven 
diffusive system s p| , as well as phase transitions far from 
equilibrium JlO| , |ll| 

The main theoretical interest in power-law scaling 
stems from the fact that the long-range properties of 
such systems are universal, i.e., they are determined by 
certain symmetry properties and do not depend on mi- 
croscopic details of the dynamics. This allows one to 
categorize phase transitions and growth phenomena into 
universality classes. Typically each universality class is 
associated with a certain set of values of the critical ex- 
ponents. Moreover, the functional form of the scaling 
function / turns out to be universal as well. 

As an example for power-law scaling, which we are go- 
ing to contrast with a different type of scaling in this 
paper, let us consider the following well-known simple 
model of molecular beam epitaxy (MBE): A particle 
beam deposits atoms onto a fiat substrate at a flux F 
(atoms per unit area per unit time). The atoms diffuse 
on the substrate with a surface diffusion constant D un- 
til they meet another adatom, in which case they form a 
stable and immobile nucleus of a two-dimensional island 
on the surface, or until they attach irreversibly to the 
edge of an already existing island. If they reach the edge 
diffusing on top of the island, they go down and attach 



to the edge with the same probability as if they arrived 
there from the lower terrace, i.e., Ehrlich-Schwoebel bar- 
riers are not taken into account. 

The observable examined in this paper is the time- 
dependent nucleation density, n, i.e., the number of nu- 
cleation events per unit area in the first layer integrated 
over time. Obviously, the nucleation density is a funda- 
mental quantity characterizing the island morphology as 
it indicates how many islands are formed. By definition, 
the nucleation density increases monotonically with time 
and saturates at a constant value when the first mono- 
layer is completed. 

In MBE the two parameters D and F can be used to 
construct a characteristic length 

£ Q = (D/F) 1 / 4 . (3) 

When the nucleation density reaches the value 1/£q, the 
rate of nucleation events decreases drastically since it be- 
comes more likely that an adatom attaches to an already 
existing island instead of meeting another adatom and 
forming a nucleus. In terms of the coverage O = Ft, i.e., 
the total number of deposited atoms per unit area, the 
time dependence of n is known to obey the scaling form 

r^ ,e)^ - 2 /(e^o)- (4) 

As shown in Ref. jl2], the scaling function / behaves as 

... ( z 3 for < z < 1 
AM«( z i/3 {oT i <<z<<i y a 2 > (5) 

where a is the lattice constant, which will be set to unity. 
Thus, MBE exhibits standard power-law scaling as de- 
scribed by Eqs. ([!]) and (||). However, in pulsed laser- 
deposition (PLD), an alternative method of growing thin 
epitaxial films, we find that the nucleation density shows 
a fundamentally different type of scaling, which in MBE 
is realized only approximately. 

In order to work out the difference most clearly, let 
us assume that the amplitudes of the asymptotic power 
laws in (^|) are equal (which is approximately the case 
in MBE), so that one can get rid of them by dividing 
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n(£o, 6) by the nucleation density at a particular cover- 
age Qq. In the following we choose ©o = 1 and define 



M(l , 8) = n(£ , e)/n(£ , 1) = f 2/3 f 2 (ef 



(6) 



The scaling function /2 is obtained from /i by re- 
placing the proportionality in (^J) by an equal-sign. If 
we furthermore suppose that the asymptotic power laws 
m (|) would remain valid right to the crossover point 
at Q£q = 1 (which is certainly not the case), then the 
scaling function would have the additional symmetry 

f 2 (z x ) = f 2 \z). (7) 

As a consequence the system would not only be invariant 
under the scale transformation (llh, 



but also under 



6 A~ 2 9, M -> A~ 2/3 M, 



to 



' o - 



(8) 



(9) 



This is a scaling transformation for the logarithms with 
all critical exponents equal to 1. According to (g) this 
would imply 



InM = (ln4) t/(ln6/ln4). 



(10) 



with a piecewise linear scaling function g(z) whose slopes 
are determined by the exponents in (^). 

While in MBE this type of scaling holds only in an ap- 
proximate sense, we are going to show that PLD indeed 
obeys this scaling form with a quadratic scaling function 
g. PLD is a growth technique in which the target ma- 
terial is ablated by a pulsed laser and then deposited 
in pulses on a substrate surface, i.e., many particles ar- 
rive simultaneously at the surface [T^ ]. Experimentally 
each pulse has a length of about a few nanoseconds and 
the time between two pulses is of the order of seconds. 
In order to investigate the scaling behavior of PLD, we 
consider a simple model which generalizes the standard 
model for MBE [jl4|jl5| . In this model the duration of a 
pulse is assumed to be zero and the transient enhance- 
ment of the mobility of freshly deposited atoms is ne- 
glected. It is controlled by three parameters, namely, the 
intensity, /, which is the density of particles deposited 
per pulse, the rate for diffusion of adatoms on the surface, 
D, and the average flux of incident particles per site, F. 
Since there is no edge diffusion, the islands grow in a frac- 
tal manner before they coalesce. A similar model with 
a finite pulse length was studied previously by Jensen, 
Niemeyer, and Combe Jl6| , p7[ . 

If the intensity / is very low, PLD and MBE dis- 
play essentially the same properties. However, if the 
intensity exceeds the average density of adatoms dur- 
ing a MBE process, we expect a crossover to a different 
type of behavior. Since this density is known to scale as 
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FIG. 1. Molecular beam epitaxy (left) compared to pulsed 
laser deposition (right) for D/F = 10 8 and I = 0.01. The 
figure shows typical configurations after deposition of 0.05 
monolayers. 



(D/F) 1 - 2 '', where 7 = 1/(4 + d f ) 
takes place at a critical intensity 



(D/F) 



-2/3 



1/6, the crossover 



(11) 



The qualitative difference between PLD and MBE for 
/ > I c is shown in Fig. [I]. As can be seen, there are much 
more nucleations at an early stage, although the effective 
flux of incoming particles is the same in both cases. 

In order to avoid the influence of the crossover at 
I ~ I c , we restrict our PLD-simulations to a particu- 
larly simple case, namely to the limit of an infinite D/F, 
meaning that all adatoms nucleate or attach to an ex- 
isting island before the next pulse arrives. In this limit 
I c = 0, and the nucleation density again depends only on 
two variables, n(I, 0). 

Performing Monte Carlo simulations we investigated 
the nucleation density for various intensities using a sys- 
tem size of 400 x 400. The measurement always takes 
place right before a new pulse is released. Our results 
are shown in Fig. |^. As can be seen, n(I, 0) increases 
with increasing intensity /. This is plausible since for a 
higher intensity more atoms arrive at the surface simul- 
taneously so that most of them can meet and form new 
nucleations before attaching to already existing islands. 

Obviously, the curves in Fig. do not display asymp- 
totic power laws, rather they are defined on finite inter- 
vals. Moreover, it is impossible to produce an ordinary 
data collapse by shifting the curves horizontally and ver- 
tically in the double-logarithmic plot. Thus, the usual 
scaling theory relying on power-law scaling fails. How- 
ever, as we are going to show below, it is possible to 
generate a data collapse by using the scaling form (|l(]) . 

The main idea is to stretch the curves in Fig. || both 
horizontally and vertically in such a way that their end- 
points collapse. To this end let us first consider the right- 
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most data point of each curve. As shown in Ref. 
the saturated nucleation density n(1, 1) scales as 

n{I,l)~P v , (12) 
where the exponent v was estimated numerically by 
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FIG. 2. The nucleation density versus time during the de- FIG. 3. The same data rescaled in such a way that all 

position of a monolayer. The dashed line has the slope 1. curves terminate at the point (1, 1). 



v ~ 0.28(1). The rightmost points can be collapsed by 
plotting the normalized nucleation density 



M(/,e) = 



n(i,e) 

n(I,l) 



(13) 



against 0, as shown in Fig. (|3|). 

Turning to the leftmost data points in Fig. (|J), which 
correspond to the measurements after the first pulse, we 
note that these points lie on a straight line with slope 1. 
This observation can be explained as follows. After de- 
position of the first pulse, the adatoms diffuse until they 
meet and nucleate. Most of them will nucleate with an- 
other adatom and form an island consisting of two atoms. 
Thus, after completion of the nucleation process, the nu- 
cleation density would be 1/2. In reality however, some 
of the adatoms form bigger islands with three or more 
particles, but to leading order these processes can be ne- 
glected. Thus, we can conclude that 



n(/,J)~/ 



(14) 



which is the vertical coordinate of the leftmost data 
points. Moreover, their horizontal coordinate of the 
points is given by = J, explaining the slope 1 of the 
dashed line in Fig. |[ In fact, the simulation data deviate 
from a linear law by less than 3 %. 

Surprisingly the power laws (|l^) and ([lj]) are valid not 
only for small I but also for I — > 1. This implies that 
the prefactors of ([l2|) and (|lj) have to be identical. This 
can be seen from Fig. where we plotted the normalized 
nucleation density M(I, 0) versus (O). In this represen- 
tation the leftmost points of the curves obey a power law 



M(I,I) 



(15) 



extrapolating to M (1, 1) = 1, which proves the equality 
of the prefactors in Eqs. (|lj) and ([hi]). 



After collapsing the rightmost data points, we collapse 
the leftmost data points by stretching the curves ver- 
tically and horizontally. This can be done by dividing 
In M and In by certain factors which are proportional 
to 1/ln/. Clearly, this manipulation does not affect the 
upper terminal point of the curves. In fact, plotting 
InM/ In/ versus ln0/ln/ (see Fig. ||), we obtain a con- 
vincing data collapse 

Assuming this data collapse to hold, the system turns 
out to be invariant under the transformation 
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x 



I" 



which implies the scaling form 



InM ~ (In/) g 



ln0 
TnT 



(16) 



(17) 



This scaling form has the same structure as the approx- 
imate scaling law (^) for MBE. Moreover, we find that 
the scaling function is well approximated by a simple 
parabola 



g[z) 



a z 



(18) 



so that InM ~ a In 2 ©/ In/. Because of the quadratic 
form the data collapse can be gener ated in several ways, 
e.g., by plotting InM versus lnO/i/hi/. 

Finally, we present a numerical puzzle. As outlined 
before, the scaling form (|l7]) for PLD should also hold for 
finite D/F provided that I > I c . Therefore, this scaling 
form can be compared with the approximate scaling form 
for MBE © at the crossover point / ~ I c , where both 
scaling concepts should 'intersect'. Because of Eqs. (||) 
and (pTp, the system is then characterized by the length 
scale £q ~ J -3 / 8 . Consequently the asympotic power- 
laws for MBE translated into the language of PLD read 
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FIG. 4. Data collapse according to the scaling form (|17|). 



Performing numerical simulations, we have demonstrated 
that such a scaling form leads to a convincing data col- 
lapse for PLD in the limit of an infinite diffusion constant. 
In contrast to ordinary scaling functions the function g 
is defined on a limited interval between two points. 

Various questions remain open. First of all, it would 
be nice to find further examples for this kind of scaling. 
In this context it would be particularly interesting to 
investigate an exactly solvable model with such a scaling 
behavior and to prove the validity of the proposed scaling 
form. Moreover, it is not yet clear to what extent the 
scaling function g is universal, i.e., independent of details 
of the dynamic rules of the model. Finally, the numerical 
coincidence of the crossover point ( pl| ) and a particular 
point of the collapsed curves is not yet understood. 

We thank the Deutsche Forschungsgemeinschaft for sup- 
port within SFB 491. 



M(i,e) 



e 3 // 2 

l/3 



for < e « i 3 / 4 

for 7 3 / 4 < 6 < 1 



(19) 



Taking the logarithm and dividing by In I we obtain 
lnM(I,e) I ; -'T^f - 2 



In/ 



olne 
d TnT 

ln6 
3 In I 



In e 
In I 



for 0< £f < f 



(20) 



In the limit I — > the crossover between the two power 
laws becomes sharper so that Eq. (^0|) converges to a 
piecewise linear curve, which is shown in Fig. | as a 
dashed line. Surprisingly, the crossover point 



ln6 c 
In/ 



= 3/4 ; 



lnA/(/,9 c ) 
In/ 



1/4 



(21) 



seems to lie precisely on the collapsed curves for PLD. If 
this were true, the prefactor a in Eq. fll8| ) would be given 
by a = 4/9. Furthermore, the scaling relation ( |i~5| ) im- 
plies that the prefactor a and the exponent v are related 
by a = 1 — 2v. Thus, we obtain the result 



(1 - a)/2 = 5/18 ~ 0.278. 



Surprisingly, this value is in perfect agreement with the 
previous numerical estimate v ~ 0.28(1) of Ref. p5| . 

To summarize, we have demonstrated that a simple 
model for pulsed laser deposition displays an unusual 
type of scaling behavior. In its most general setup, this 
kind of scaling behavior is observed in systems which are 
invariant under the transformation 



M -> M 



A" 



for arbitrary A, where a and /3\, . . . , /3„ are certain expo- 
nents. The corresponding scaling form reads 



lnM~ (lnzi) Q 5 



lnz 2 
(In 21^ 



In z n 



(lnzi)/ 3 - 
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